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Alongside the effort underway to build quantum computers, it is important to better understand 
which classes of problems they will find easy and which others even they will find intractable. We 
study random ensembles of the QMAi -complete quantum satisfiability (QSAT) problem introduced 
by Bravyi [l]. QSAT appropriately generalizes the NP-complete classical satisfiability (SAT) prob- 
lem. We show that, as the density of clauses/projectors is varied, the ensembles exhibit quantum 
phase transitions between phases that are satisfiable and unsatisfiable. Remarkably, almost all in- 
stances of QSAT for any hypergraph exhibit the same dimension of the satisfying manifold. This 
establishes the QSAT decision problem as equivalent to a, potentially new, graph theoretic problem 
and that the hardest typical instances are likely to be localized in a bounded range of clause density. 
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I. INTRODUCTION 

The potential power of quantum computers motivates 
the intense effort in progress to understand and, even- 
tually, build them. Much interest has, naturally, been 
focused on algorithms that outperform their classical 
counterparts. However the development of a complex- 
ity theory for quantum computers suggests that we al- 
ready know problems, those shown to be QMA-complete, 
whose worst case solutions will take even quantum com- 
puters a time exponential in problem size. These appro- 
priately generalize the class of NP-complete problems- 



those which are easy to check but (believed to be) hard 
to solve on classical computers-to quantum computers 

mm- 

The classic technique of complexity theory is assign- 
ing guilt by association, i.e. by exhibiting the polyno- 
mial time equivalence of a new problem to a particular 
problem believed not to be amenable to efficient solution. 
For NP-complete problems this leads to the celebrated 
Cook-Levin theorem which shows that the satisfiability 
problem with clauses in three Boolean variables (3-SAT) 
encapsulates the difficulty of the entire class [4 . While 
this is a powerful and rigorous approach it has two lim- 
itations. It does not directly tell us why a problem has 
hard instances and it does not tell us what general fea- 
tures they possess. 

Over the past decade joint work by computer scien- 
tists and statistical physicists has produced an interest- 
ing set of insights into these two lacunae for 3-SAT (and 
related constraint satisfaction problems). These insights 
have come from the study of instances chosen at ran- 
dom from ensembles where the density of clauses a acts 
as a control parameter. One can think of this as repre- 
senting the study of typical, but not necessarily worst, 
cases with a specified density of clauses/constraints. Us- 
ing techniques developed for the study of random systems 
in physics, it has been shown that the ensembles exhibit a 
set of phase transitions between a trivial satisfiable (SAT) 
phase at small a and an unsatisfiable (UNSAT) phase at 
large a (see e.g. [3 El 13 [H [H])- These transitions mark 
sharp discontinuities in the organization of SAT assign- 
ments in configuration space as well a vanishing of SAT 
assignments altogether. This information has provided 
a heuristic understanding of why known algorithms fail 
on random SAT when the solution space is sufficiently 
complex and in doing so has localized the most difficult 
members of these ensembles to a bounded range of a. 

In this work we begin an analogous program for quan- 
tum computation with the intention of gaining insight 
into the difficulty of solving QMA-complete problems. 
Specifically, we introduce and study random ensembles of 
instances of the quantum satisfiability (fc-QSAT) prob- 
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lem formulated by Bravyi [T]. Like classical 2-SAT, 2- 
QSAT is efficiently solvable (i.e. it is in P), while for 
/c > 4, fc-QSAT is QMAi-complete [TO]. 

After laying out some important definitions and back- 
ground on our problem in Sec.|ll] we attack the ("easy") 
k — 2 problem in Sec. Here, we derive the complete 
phase diagram using transfer matrix techniques intro- 
duced by Bravyi In the process, we discover that 
the classical geometry of the 2-QSAT interaction graph 
determines the generic satisfiability of random instances 
without reference to the random quantum Hamiltonian 
imposed upon the graph. We further show that our ran- 
dom ensemble asymptotically satisfies a technically im- 
portant constraint known as the promise gap with prob- 
ability exponentially close to 1. 

In Sec.|TV]we move on to establish rigorous bounds on 
the existence of SAT and UNSAT phases for the ("hard" ) 
k > 3 problem. Although the proof is quite different 
from the analysis for fc = 2, we then directly recover the 
unexpected deduction that the satisfiability of a generic 
instance of fc-QSAT reduces to a classical graph prop- 
erty. This allows us to more tightly bound the UNSAT 
phase transition. We comment briefly on the straight- 
forward generalizations to some related ensembles and 
on the satisfaction of the promise gap. Finally, in Sec. |V] 
we lay out the goals for the continued study of random 
quantum problems and algorithms. 



II. DEFINITION OF RANDOM QSAT 

We consider a set of N qubits. We first randomly 
choose a collection of fc-tuples, {Im,rn — 1 . . . M} by 
independently including each of the (^) possible tuples 

with probability p — a/(^). This defines an Erdos hy- 
pergraph with aN expected edges; we exhibit simple ex- 
amples of these for fc = 2, 3 in Fig. [T| In classical fc-SAT, 
the next step is to generate an instance of the problem 
by further randomly assigning a Boolean k-clause to each 
k-tuple of the hypergraph. In key contrast to the clas- 
sical case, where the Boolean variables and clauses take 
on discrete values ~ true or false - their quantum gen- 
eralizations are continuous: the states of a qubit live in 
Hilbert space, which allows for linear combinations of 
|0) ("false") and |1) ("true"). Thinking of a Boolean 
clause as forbidding one out of 2^ configurations leads to 
its quantum generalization as a projector IlJ, = 
which penalizes any overlap of a state {ip) of the fc qubits 
in set / with a state |0) in their 2'^ dimensional Hilbert 
space. In order to generate an instance of fc-QSAT the 
states 10) (of unit norm) will be picked randomly in this 
Hilbert space. 

The sum of these projectors defines a positive semidefi- 

M 

nite Hamiltonian H = ^ H^"* and the decision problem 

771 — 1 

for a given instance is, essentially, to ask if there exists 
a state that simultaneously satisfies all of the projectors. 



i.e. to determine whether H has ground state energy, 
Eo, exactly zero. The quahfier "essentially" is needed 
because Eq is a continuous variable and in order that the 
problem be checkable by a quantum verifier (and there- 
fore in QMAi), H must be accompanied by a promise 
that Eo is either exactly zero or exceeds an e ~ 0{N^°-) 
where a is a constant. The scaling function e is also 
known as the promise gap. For our random ensemble we 
wish to compute the statistics of this decision problem 
as a function of a. Specifically we would like to know if 
there are phase transitions, starting with a SAT-UNSAT 
transition, in the satisfying manifold as a is varied. Ad- 
ditionally, we would like to check that the statistics in 
the large A'^ limit are dominated by instances that auto- 
matically satisfy the promise gap. 

We now note a few key properties of the hypergraph 
ensemble: At small clause density a, the size distribu- 
tion of connected components ("clusters") is exponential. 
Moreover, almost all of these clusters are treelike with a 
few, 0{N'^), containing a single closed loop (see Fig. [ij; 
the number of clusters with several closed loops vanishes 
as —> oo. Above a critical value agc(fc) = -jTQ^zr^j ^ 
giant component emerges on which a finite fraction of 
the qubits reside. Unlike the finite clusters, the giant 
component may contain a non- vanishing (hyper)core, de- 
fined as the subgraph remaining after recursively strip- 
ping away leaf nodes {i.e. nodes of degree 1). For 
fc = 2, an extensive hypercore emerges continuously at 
Q!hc(2) = agc(2) = |; for fc > 3, the hypercore ap- 
pears abruptly with a finite fraction of the nodes at 
ahc{k) > agcik) [mmi]. 

III. 2-QSAT 

Let us first consider our questions in the specific con- 
text of 2-QSAT (phase diagram in Fig. |2]). While this 
is a classically easy (P) problem, the random ensemble 
has much of the structure we also find for the harder 
fc > 2 case. An instance H of QSAT is satisfiable if its 
kernel ker(i?) has nonzero dimension; in the following, 
we will find the kernel of H explicitly by transforming 
the random projector problem almost surely into a sim- 
ple Heisenberg ferromagnet, whose ground state space is 
well known. 

A central tool in this analysis is Bravyi 's transfer ma- 
trix T^^ which, given a vector 1-0^) of the state of qubit 

i yields \ipj) = T^^'\'^pi) for qubit j, such that the prod- 
uct state \tpij) = (8) satisfies the projector onto 0: 
(i/'ij |H^ IV'ij) = 0. Concretely, if the projector penalizes a 
joint state of both qubits given by the complex vector = 

(000,001,010,011), we have = e0t = (+^2^ +^^1). 
Here e is the standard antisymmetric matrix (Levi-Civita 
symbol) in two dimensions. We note that the transfer 
matrix T*-' for any given link is almost surely invertible 
so this construction finds satisfying product states for 
any input \ipi). 
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FIG. 1: Examples of random graphs and hypergraphs for (a) 2-SAT and (b) 3-SAT, respectively, (a) The clusters, clockwise 
from bottom left, are chain, tree, clusters with one and two closed loops ("figure eight"). The short closed loops, as well as 
the planarity of the graphs, are not representative of the large- limit. The two-core of the graphs (shaded blue) is obtained 
by removing the (unshaded, green) dangling tree structures, (b) Each square represents a hyperedge connected to 3 nodes. 
Clockwise from top left are a tree, a hypercore and a hypergraph with simple loops (shaded red) but vanishing hypercore. 



A. Trees are SAT 

Consider the clusters that enter the 2-QSAT graph en- 
semble. Of these, a tree comprising n qubits and n — 1 

n 

edges has a satisfying product state \^) = where 

is obtained from an arbitrary reference qubit i — 1 
by repeated application of the T's along the (unique) 
path joining i with j. In fact, the satisfying subspace 
is, almost always, n -I- 1 dimensional. To show this, we 
will map the random projector problem directly onto the 
ferromagnetic Hamiltonian on the same tree. In more 
mathematical terms, we construct a non-unitary action 
of the permutation group on the qubit Hilbert space that 
leaves the zero energy space ker(iJ) invariant. Thus, 
the ground state space is precisely the totally symmetric 
space Sym"C^, which has dimension n-\- \. 

We define a particular, non- orthogonal product basis 
for the Hilbert space H = H° «) «)•••«) H""^ First, 
choose any two linearly independent unit vectors [t*^) and 
ll'^) as a basis for TL^ . Second, use the transfer matrix 
to transfer this basis of to a normalized basis 



\v) 



(1) 



of W for each of the neighbors j of 0. Finally, recursively 
traverse the tree to produce a pair of vectors |t*) , 
for each site i in the tree. This procedures produces a 
choice of basis for each of the individual qubit Hilbert 
spaces which we use to define a product basis we call 



the transfer basis. Although we use the notation of spin 
up and down, we emphasize that the states need not be 
orthogonal. 

Finally, we show that the ground state space of the 
Hamiltonian H is precisely the space of totally symmet- 
ric states defined with respect to the transfer basis. We 
consider the constraint that a generic vector e 7i is 
annihilated by H"^ by factoring 



+ (|fi^)-|i"T^)) 



(2) 



The first three terms are identically annihilated by the 
projector H"^. This follows trivially from the construc- 
tion of the transfer basis for the first two terms while for 
the third term, keeping track of indices carefully (sum- 
mation over repeated indices is understood): 

n°H|TT) + |ir)) = |r)((0°^*U.a. CM 

+ )ao,ai iagtai ) 







(3) 



where we have exploited the antisymmetry of e. Thus, 
|u4) must be zero and ^! clearly lies in the symmetric 
eigenspace for swaps (01). A similar argument holds for 
each of the swaps {ij) on links of the tree B. Since B 
is connected, these swaps generate the full permutation 
group Sn and \E' must be in the completely symmetric 
subspace of its action. In particular, this means that 
the ground state space is isomorphic to Sym"C^ which is 
n + I dimensional. 
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B. Loops 

To understand hypergraphs with loops, we first look 
for satisfying product states. For a cluster G compris- 
ing TiL independent closed loops {i.e. n qubits with 
n + TiL — 1 edges, see Fig. [T|, a product state is inter- 
nally consistent only if the product of the T's around 
each closed loop returns the state it started with. For 
a graph with a single closed loop, Oi, this imposes 

(rijfceOi T^^^ \t\jij = AilV'i), which in general allows two 
solutions, so that these graphs are SAT. By contrast, if 
a site i is part of a second independent closed loop, O2, 

the additional demand (^njfce02 -^"''°) l''^^) ^ '^2|V'j) al- 
ready yields an overconstraint. Thus graphs with more 
than one closed loop lack satisfying product states with 
probability 1. 

In fact, we can generalize the above reasoning to show 
that it holds even for the entangled (non-product) states 
in the problem as follows. We choose a spanning tree 
on G and starting qubit to define a transfer basis - 
the ground state space is necessarily a subspace of the 
associated symmetric space. In particular, if G contains 
a single closed loop Oi, we take qubit on the loop. 
In this case, we choose the starting basis on qubit to 
satisfy the product state consistency condition for 0\ 



n 



(4) 



i.e. we take the eigenbasis of the loop transfer matrix, 
which we denote t"ii" with eigenvalues A|,A|. These 
two states transfer to two linearly independent product 
ground states for the loop - all up and all down. A short 
calculation verifies that in fact there are no more: write 
a generic ^ G ker(ff) C Sym^V'^ 



as 



i*> = \r'f)\vi) + \r\n\v2) 
+ {\r'n + \rr))\v,), 



(5) 



where —1 is the qubit at the end of the loop Oi and Vi 
are vectors on the n — 2 other qubits. The projection 
lli-ifi) annihilates the first two terms by choice of the 
basis. However: 



n 



(-1,0) 



0(-l.O>^(|-lT ^(-l,0>t ^0 

+h l"^ ^ f) 1-3) 
'^<-^'°>)(AT-Ai)(r6f)|.3) (6) 



Since A| Aj^ w.p. 1, V3 must be zero and by symmetry, 
ip must be in the span of the all up and all down states. 

An analogous calculation verifies that any additional 
closed loops introduce projectors that are violated on 




/ 0{a — ac)^ 



FIG. 2: Phase diagram of 2-QSAT. The quantum SAT- 
UNSAT transition a, = | coincides with the emergence of 
a giant component in the random graph, which lies at half 
the classical 2-SAT transition Uc — 1 12:. The solid (dashed) 
line marks an asymptotic upper bound on the quantum (clas- 
sical) ground state energy density £o- 



this two-dimensional subspace and thus graphs with more 
than one loop are unsatisfiable. 

It is worthwhile to consider how the usual Heisenberg 
ferromagnet fits into the above calculations. In this case, 
the transfer matrices are all the identity and thus all of 
the loop constraint conditions are identically satisfied. 
From the point of view of Eq. [6] the ferromagnet has 



At 



A, 



1 and there is no reduction in the ground 



state degeneracy due to the introduction of loop-closing 
ferromagnetic bonds. 



C. Phase diagram 

In light of the above analysis, the existence of a 
SAT/UNSAT phase reduces to the presence of multiply 
connected components in the ensemble of 2-SAT graphs 
(rather than the combined ensemble of graphs and pro- 
jectors). For a < age = 5, the number of such clusters 
vanishes in the limit of large N, so any instance is SAT 
with probability that goes to 1 as ^ 00. At age = 1/2, 
closed loops proliferate as a giant component appears and 
thereafter all instances are UNSAT with probability 1. A 
straightforward upper bound on the energy in the UN- 
SAT phase, E < 0{{a - agc)'^N), follows from the fact 
that the fraction of sites in the core of the giant compo- 
nent grows as (a — ccgc)^. 

Physical intuition suggests that the ground state en- 
ergy should be extensive above the transition and thus 
likely saturate the given upper bound. In the next sec- 
tion, we adduce strong evidence for this by exhibiting a 
minimal lower bound on the energy E > iV^"*^ for any 
e > which holds except for exponentially rare instances. 
This follows from the twin claims that (i) the energy of a 
figure eight comprising n sites decays only polynomially 
with n, and (ii) that the number of disjoint figure eights 
in the random graph grows nearly linearly with N. Ob- 
serve that this kind of lower bound is exactly what we 
need to establish that our ensemble keeps the promise 
that either E = or E > \/N°' with probability ex- 
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ponentially close to 1. While this demonstration is not 
strictly needed for 2-QSAT (since it is in P) , it is sugges- 
tive of what we might expect for the k > 3 cases. 



D. The promise gap for k — 2: counting figure 
eights 

We expect on physical grounds that the UNSAT phase 
of fc-QSAT has extensive ground state energy with rela- 
tively vanishing fluctuations for any k. In this case, the 
promise that E > 0(7V~°) fails to be satisfied only with 
exponentially small probability by Chebyshev's inequal- 
ity. More generally, so long as the average ground state 
energy is bounded below by a polynomially small scale 
E > 0{N~^) with relatively vanishing fluctuations, the 
promise will be violated with only exponentially small 
probability for a > b. 

Placing rigorous lower bounds on the expected quan- 
tum mechanical ground state energy is generally diffi- 
cult, but at least for fc = 2, we argue as follows to flnd 
a nearly extensive bound scaling as N^~'^ for any e > 0. 
As the Hamiltonian for QSAT is a sum of nonnegative 
terms, we can bound the ground state energy from below 
by considering manageable subgraphs and ignoring the 
contribution from other terms. In particular, the aver- 
age ground state energy of a figure eight graph, that is 
a loop of length L with one additional crossing edge, is 
polynomially bounded below by 0{L~^). This already 
gives a polynomial lower bound on the expected energy 
simply by allowing L to scale as N and knowing that we 
will find at least one such subgraph in the giant compo- 
nent with probability exponentially close to 1. We will 
do better by finding a large number of disjoint such figure 
eight graphs. 

The expected number of subgraphs A in the random 
graph Gjv.p is given by the following formula: 

iV! 



E#(A c G) = 



,(A) 



(7) 



{N-\A\)\\Aut{A)\ 

where \A\ is the number of vertices in A, e{A) is the 
number of edges in A and Aut{A) is the group of auto- 
morphisms of A {\Aut{A) \ is the cardinality of this set). 
This formula simply counts the number of ways of find- 
ing permutations of \A\ nodes in G and connecting them 
up into an A subgraph. For the fixed clause density en- 
semble, we take p = 2a/N. 

A figure eight graph is uniquely specified by giving its 
size L (we take L even) and the distance d = 2 . . . L/2 
between the two nodes that are connected by the crossing 
link. We let A be the disjoint union of K figure eight 
graphs with L nodes each and cross bars at separation 
L/2 - 1. Such a graph has \Aut{A)\ = K\2^ from the 
Kl permutations of the disjoint subgraphs and the two- 
fold symmetry of each figure eight. Thus, the expected 
number of K-iold disjoint figure eights of size L is 

]Sf\pKiL+l) 



S^KLl log 2a 



We now allow if, L to scale with TV such that KL <C N 
and use Stirling's formula to find the asymptotic entropy 

KL \ 

— j +K (log a + 1 - log i-T - log iV) . 

(9) 

This entropy is positive and growing with TV so long as 
L log A^, KL <C N and a > = 1/2. In this regime, 
assuming that the fluctuations in S are relatively small, 
we flnd that there are K ^ N^~'^ disjoint flgure eights of 
size L ^ N'^/'^ with exponentially high probability. These 
lead to a nearly extensive lower bound on the expected 
ground state energy. 



IV. fc-QSAT 

We now turn to fc-QSAT with k > 2. As noted above, 
for sufficiently small a, the important hypergraphs have 
vanishing cores. In this regime, it is possible to construct 
the hypergraph by adding in edges one by one where 
each additional edge brings with it at least one new leaf 
node. Thus by a generalization of the transfer matrix ar- 
guments for A; = 2, it is possible to construct a satisfying 
product state on the full hypergraph; the details are in 



E# 



{N - KLy.K\2^ 



(8) 



Sec. IV A Deducing the existence of the UNSAT phase 
is not as straightforward however for, unlike in the k = 2 
problem, we do not know a set of unsatisflable graphs 
which are present with flnite probability as A'^ — > oo; in- 
deed these are not known for classical SAT either. Thus, 
in Sec. |IVB| we produce an indirect proof that the di- 
mension of the kernel of H vanishes for sufficiently large 
a. 

Recall that for 2-QSAT the dimension of the satisfying 
manifold ended up depending, with probability 1, only 
on the topology of the graph and not on the choices of 
the projectors. This is in fact a very general result, which 
we prove directly in Sec. |IVC| For any hypergraph for k- 
QSAT (regardless of its likelihood) the dimension of the 
satisfying manifold is independent of the choice of pro- 
jectors with probability 1. Moreover, non-generic choices 
of projectors result in larger satisfying manifolds. Implic- 
itly, this means that satisflability for the random ensem- 
ble is a graph theoretic property and raises the extremely 
interesting possibility that random /c-QSAT can be for- 
mulated without reference to quantum Hamiltonians at 
all. 

Finally, we comment on the satisfaction of the promise 
gap in Sec. |IVD and on the generalization of QSAT to 
higher rank projectors in Sec. |IVE[ 



A. Existence of SAT phase 

We prove the existence of a satisflable phase at low 
clause density for rank 1 fc-QSAT by constructing prod- 
uct states of zero energy for arbitrary hypergraphs con- 
taining a hypercore with a satisfying product state. Since 
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the core of a random graph vanishes for a < ahc(fc), this 
proves the existence of a satisfiable phase below the emer- 
gence of a hypercore. 

Lemma 1. Suppose that H is an instance of QSAT on N 
quhits with a satisfying product state \^) . Let H' = 7? -1-11 
be an instance of QSAT on N' > N quhits where 11 is a 
QSAT projector touching at least one quhit not among 
the original N . Then H' has a satisfying product state 

Proof. Assume first that the hyperedge 11 adjoins only 
one dangling qubit (untouched by H), say qubit 0. Let 
the components of \^) be The key observation 

here is that the projector H with fc— 1 of its qubits fixed in 
a product state of N qubits uniquely specifies the state 
of the dangling qubit 0. If 0.^^...^^, ^ is the vector onto 
which n projects, then the state of qubit must be per- 
pendicular to '?!>*o...,,_i^'ii...ifc_i. 

Thus, contraction with e*o*° defines a fc — 1 to 1 transfer 
matrix T^°...j^^ ^ generalizing the Bravyi transfer matrix 
of 2-QSAT. The new product state given by l^*') = 
\^/) where i>°^ = ^/"...j^^i^Ji -ifc-i satisfies n. Moreover, 
since H does not act on qubit 0, H \^') = |*) = 

0. Thus H' = H + U is satisfied by the product state l^*') 
on + 1 qubits. 

Finally, if H adjoins more than one dangling qubit, 
simply fix all but one of these to an arbitrary state and 
then apply the above transfer matrix procedure. □ 

Theorem 2. Suppose that H is an instance of random 
QSAT and let He be the restriction to its hypercore. If 
there is a satisfying product state on He, then there is 
almost surely a satisfying product state on H. 

Proof Let n,, i = 1 • • • A/ (G) - M (C) be the sequence of 
projectors removed in the process of stripping leaf hyper- 
edges from the hypergraph of H to produce the hyper- 
core. At each stage of this process, the removed projector 
Hi has at least one dangling qubit. Thus, if we iteratively 
reconstruct H from He by adding back the Hi in reverse 
order, we can apply the lemma at each stage to lift the 
product state on He to a product state on H. □ 

B. The UNSAT phase exists 

For a given random hypergraph, let us construct H via 

m J 

the sequence H„i = J2 LI^ . , i-e. we add the projectors 

one at a time in some order. Clearly H ~ H^. Let 
Dm be the dimension of the satisfying subspace for Hm', 
evidently Dq — 2^ . At each step, if the next added 
projector involves a set of qubits that were not previously 
acted upon, then Dm+i — Dm{^ — ^) as we may simply 
implement the projection by reducing the size of a basis 
that can be factorized between the target qubits and all 
others. It is intuitively plausible that this is the best we 



can do — in cases where the target qubits are entangled 
with the rest of the system we should expect to lose even 
more states, i.e. 

in general. A proof that this bound holds with proba- 
bility 1 for projectors randomly chosen according to the 
uniform Haax measure is contained in Appendix|X] From 
this we conclude that 

Thus, for a > -l/log2(l- ^) ~ 2*" the problem is 
asymptotically almost always UNSAT. 

C. Geometrization theorem 

Geometrization Theorem. Given an instance H of 
random k-QSAT over a hypergraph G, the degeneracy of 
zero energy states dim(ker(if)) takes a particular value D 
with probability 1 with respect to the choice of projectors 
on the edges of the hypergraph G. 

Proof. For a fixed hypergraph G with M edges, H = 

^4> = Y^iii = Y^fLi {(pil is a matrix valued func- 
tion of the components of the set of M vectors \(j)i). 
In particular, its entries are polynomials in those com- 
ponents. Choose |</)) such that H has maximal rank R. 
Then there exists a,n R x R submatrix of H such that 
det(H\f(xR.) is nonzero. But this submatrix determinant 
is a polynomial in the components of |0) and therefore 
is only zero on a submanifold of the |0) of codimension 
at least 1. Hence, with probability 1, H has rank R and 
the degeneracy dim(ker(_ff)) = 2^ — i?. □ 

The theorem holds for general rank r problems as well 
by a simple modification of the argument to allow extra 
(j)'s to be associated to each edge. 

A nice corollary of this result is a more stringent upper 
bound on the size of the SAT phase. Consider any assign- 
ment of classical clauses on a given hypergraph: we can 
also think of this as a special instance of fc-QSAT where 
the projectors are all diagonal in the computational basis. 
As this is a non-generic choice of projectors, the dimen- 
sion of its satisfying manifold is an upper bound on the 
dimension for generic choices. We conclude then that 
the classical UNSAT threshold is an upper bound on the 
quantum threshold. Indeed, if we can identify the most 
frustrated assignment of classical clauses on a given hy- 
pergraph, i.e. the assignment that minimizes the number 
of satisfying assignments, we could derive an even tighter 
bound. 

Corollary 3. The zero state degeneracy of the random 
quantum problem is bounded above (w.p.l) by the number 
of satisfying assignments of the most constrained classical 
k-SAT problem on the same graph. 
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SAT UNSAT 
0.17 0.81 4.26 7.49 

^ ^ ^ ^ 

Ogc ahc Qcav Qwb a 

FIG. 3: Phase diagram of 3-QSAT. For a < ahc « 0.81 
[TT] . the problem is rigorously SAT. The UNSAT transition is 
bounded below by the weak bound Qf„b ~ — l/log(l — 1/2*^), 
but also (shghtly less rigorously) by the classical cavity tran- 
sition Qcav [7] • The giant component emerges squarely in the 
SAT phase at Qgc = ~ 0.17. 



We have encountered many numerical examples in 
which the quantum problem has fewer ground states than 
the most frustrated classical problem on the same hyper- 
graph. Thus, the bound of corollary [3] is not tight. 

In the k = 2 case, corollary [3] gives us another way 
to show the critical point corresponds to the emergence 
of a two-core: once a two-core exists, it will contain a 
giant loop with two crossing bonds. This graph can be 
made unsatisfiable by setting each of the clauses around 
the loop to disallow the state 01 and the crossing bonds 
to disallow 00 and 11 respectively. Thus, the quantum 
problem is UNSAT w.p. 1. Ideally, we could extend this 
geometric characterization to the k > 2 problem, but we 
have thus far failed to identify a family of unsatisfiable 
hypergraphs that appear in the random hypergraph with 
non- vanishing probability. 



D. Satisfying the promise 

As in the case of 2-QSAT, physical experience sug- 
gests that the ground state energy for fc-QSAT should be 
extensive with small fluctuations above the satisfiability 
transition. If this is true, the promise is satisfied with 
probability exponentially close to 1, as in the k — 2 case. 
Unfortunately, we do not know a set of unsatisfiable sub- 
graphs analogous to the figure eight's that might be used 
to provide a more rigorous bound. 



E. QSAT at rank r > 2 

We now briefiy consider the extension of fc-QSAT to 
{k, r)-QSAT, in which the projectors 11^ have rank r, i.e. 
they penalize a uniformly chosen r-dimensional subspace 
of the 2*^ dimensional k qubit Hilbert space. The main 
results regarding fc-QSAT generalize naturally: satisfia- 
bility is almost surely only dependent on the underlying 
graph and the weak bound on the ground state degener- 
acy becomes -D^*^ < 2^ (l — , implying a bound 
ac'^^ < — l/log2 (l — ^). However, there need not be a 
SAT phase at all: if r = 2*^, the projectors 11^ are each 
the identity and the ground state energy is aN > for 
any positive a. More generally, there is some critical rank 
Tc above which the SAT phase disappears. 



We bound Tc above by 2*^/2 by exhibiting an unsat- 
isfiable subgraph of the random hypergraph that arises 
asymptotically almost surely (N — > oo) in the random 
graph ensemble even at small a. Consider a chain with 
two clauses and one shared qubit. Classically, this corre- 
sponds to a two-clause problem on 2fc — 1 bits where we 
allow each clause to forbid r = 2^ 12 configurations. Now 
let the first clause forbid all configurations in which the 
shared bit is and the second clause all configurations 
in which it is 1. This classical problem is clearly un- 
satisfiable and therefore so is any quantum problem on 
this subgraph w.'p. 1. Indeed, since there are extensively 
many such small chain components in the hypergraph, 
each with an independently chosen 0{N°) ground state 
energy, this provides an extensive lower bound on the 
total ground state energy. 

The bound < 2*^/2 is not tight however. For ex- 
ample, for fc = 3, an open chain of length four with 
rank 3 < 2^/2 = 4 is quantum mechanically unsatisfi- 
able w.p. 1, as can be checked numerically, and therefore 
so is the ensemble for (fc,r) = (3,3). However, there is 
no classically unsatisfiable problem on this chain and it 
is harder to construct a rigorous bound that scales with 
fc using this starting point. 



V. CONCLUSIONS AND OUTLOOK 

As in the classical work that was our inspiration, we ex- 
pect that the phase structure that we have identified also 
shows up in the actual performance of potential quan- 
tum algorithms. Specifically, the classical results and the 
general intuition from the study of quantum phase tran- 
sitions [13l [14] both strongly suggest that the decision 
problem will be hardest near the SAT-UNSAT bound- 
ary and easy away from such a bounded range. This 
expectation is given further force by the graph theoretic 
formulation of typical instances whereby we may again 
expect that graphs with an intermediate clause density 
are the ones where it might prove hardest to decipher 
the relevant property. As a first step in testing this ex- 
pectation, we have begun an investigation of the simplest 
adiabatic algorithm for this problem |15j and the relevant 
excitation spectra. 

Looking ahead, the immediate challenge is to explic- 
itly identify the graph theoretic property that character- 
izes satisfiable instances of random QSAT without ref- 
erence to the quantum mechanical problem. If this can 
be accomplished it will shed light on the difficulty of de- 
ciding generic instances of random QSAT and provide 
a surprising connection to classical complexity classes. 
Separately, we intend to investigate whether there are 
analogs of the clustering transitions of classical fc-SAT in 
the quantum problem and whether generalizations of the 
cavity method to quantum problems [16l [ITl [18] can help 
locate them. Finally, we would like to translate this im- 
proved understanding of this ensemble of problems into 
new algorithms for solving them on the lines of the be- 
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APPENDIX A: WEAK UNSAT BOUND 

Weak UNSAT Bound. Given an instance of random 

{k,r)-QSAT with Hamiltonian H over a hypergraph G 
with N qubits and M clauses, the degeneracy of zero 



energy states (i.e. the number of satisfying assigments) 
D = dim(ker(i7)) is bounded above by 

with probability 1. 

Proof. For a given random hypergraph G, we construct 
H = Hm via the sequence Hm = X^JLi where m = 
1 . . . M. That is, we add projectors one at a time in some 
fixed order. Let be the degeneracy of zero energy 
states after m projectors have been added. Clearly, 

Do = 2^. (A2) 

The result will follow by induction after we show that 
each additional projector reduces the degeneracy by at 
least a factor 1 — ^ w.p. 1. That is, 

Dm+l < -Dm (l - ^) . (A3) 

We consider adding a projector 11 = ^^0=1 \^") 
where = Sa/s, to the Hamiltonian H to construct 

H'. By appropriate reordering, we can always assume 
that n acts on the first k qubits of the full N qubit Hilbert 
space. We choose a basis for the zero energy subspace 
ker{H) prior to the addition of 11 as follows 

L 

\d) = T.H)\l') (A4) 

i=l 

where d € {1,. ..,£>} labels the basis states, |a^) G <C^'° 
are vectors in the first k qubit factor of the Hilbert space 
and are L linearly independent vectors on the re- 
maining {N — k) qubit fac;tor. 

The kernel oi H' = H + 11 is the subspace of the kernel 
of H which is annihilated by H. We write a generic vector 
of kev{H) as: 

D 

|V') = ^V'd|d). (A5) 
Annihilation by H leads to the condition 

= H = in E (-^"i^rf) (A6) 

i=l a=l d 



which by linear independence of the \(lf) |P) requires, for 

any a, i: 

D 

= ^<</)«K)V<i. (A7) 

That is, must he in the kernel of the rL x D matrix 
^(ia)d = ('/'"kd)- We now claim that has rank R 

at least with probability 1, which will prove our 

inductive step. This follows from two observations: 

1. The rank i? of A is a bounded random variable that 
takes on its maximal value over the choice of |^") 
with probability 1. A can be viewed as a matrix 
of monomials in the components of | Choose 
an orthonormal frame \(f)°') maximizing the rank R 
of A\ with this choice, there exists some R x R 
submatrix of A such that det(^|/jxfi) is nonzero. 
But this submatrix determinant is a homogenous 
polynomial in the components of |^") and therefore 
is only zero on a submanifold of codimension at 
least 1 in the complex Grassmannian Gr(r, 2'^) (i.e. 
the space of choices of rank r projectors 11) . Hence, 
almost every matrix A will have maximal rank R. 

2. We now need only exhibit one set of \(f)°')'s such 
that A has rank R > . This is easy: at least 
one i?-tuple of the standard basis vectors will pro- 
vide such an A. Consider the matrix of vectors |a^) 
and let ri, r2, • • • , be the ranks of each of the 2^ 
component matrices (ze. the matrices A obtained 
by using the standard basis elements for |^)). The 
Tj are the number of linearly independent rows in 
each of these matrices. Now concatenate each of 
these component matrices vertically into a giant 
2'^L X D matrix. The row rank of this matrix can 
be no larger than ri + r2 + • • • + r2fc by construc- 
tion, but there must be a full D linearly indepen- 
dent columns: any linear relation on the columns 
of the giant matrix corresponds to a linear relation 
Ylid''^d l^d) ~ ^ which Ufts trivially to a linear re- 
lation among the basis vectors \d). Hence, 

r\+r2-\ h r2fe >D. 

From this relation, there must exist some collection 

n- of size r such that ri^-\ h r^^ > rD/2''. This 

collection of basis vectors provides the frame we 
desire. 



□ 



